Modelling the onset of oxide formation 
on metal surfaces from first principles 
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The formation of ultrathin oxide layers on metal surfaces is a non-thermally-activated 
c^ ' process which takes place spontaneously at very low temperatures within nanoseconds. 

'^ ■ This paper reports mechanistic details of the initial oxidation of bare metal surfaces, 

O ! in particular Al(lll) and TiN(OOl), as obtained by means of first-principles molecu- 



lar dynamics modelling within the Density-Functional Theory. It is shown that the 



^ ■ reactions of bare metal surfaces with O2 molecules take place according to a "hot- 

o: 

en . atom" dissociative mechanism which is triggered by the filling of the a* antibonding 

cn I molecular orbital and is characterised by a sudden release of a large amount of kinetic 

t^ ■ 

Q ' energy. This released energy provides a driving force for metal/oxygen place-exchange 

00 : 

O ' processes which are responsible for the onset of oxide formation at virtually K and at 

-^ , oxygen coverages well below 1 monolayer (ML). Further simulations of the oxidation 

^ ' reactions reveal that a disordered ultrathin oxide forms on Al(lll), whereas a rather 
^ '. 

ordered structure develops on TiN(OOl) following a selective oxidation process which 

leaves clusters of Ti vacancies in the TiN lattice underneath the oxide layer. 



1. Introduction 

Metal surfaces in contact with the atmosphere become quiclky covered by oxide layers 
of thicknesses which vary from few atomic layers (less than 1 nm) to several microns, 
depending on the reactivity of the metal and the history of the sample [1]. Even the 
noblest metals, such as gold, may chemically interact with oxygen species [2] at their 
surfaces. The oxide skin governs the chemical and physical interactions of the metal 
with the outer environment, thus a precise knowledge of its structure and composition 
is fundamental to a number of technologically important problems. These comprise for 
instance the behaviour of metallic biomedical implants within the chemically aggressive 
physiological environment [3], the response of electrodes or sensors in direct contact 
with liquid solutions [4, 5] , the performances of metallic contact structures in electronic 
devices, the wear and friction properties of metal surfaces [6] , the rate of sub-critical 
crack propagation due to stress-corrosion phenomena [7], and many others [8]. 

Traditionally, engineers have been much concerned with oxidation phenomena at 
high temperatures, since they result in a continuous consumption of metal and therefore 
inevitably limit the life of machine components such as, for instance, turbine blades [1]. 
For this reason, great effort has been expended in the last century to develop a basic 
understanding of the kinetic laws governing high-temperature oxidation, and several 
theories have been formulated which explain and predict the oxidation behaviour of 
many metals and metal alloys [9-11]. However, as noted in [1], many of the existing 
rate laws seem to be satisfactory only over a limited oxide thickness range, and often 
are not transferable to low temperature or low pressure regimes. 

In particular, detailed knowledge of processes of oxide formation at room tempera- 
ture or below is scarce compared with the amount of information available in the field 
of high-temperature oxidation. The available kinetic models rely on assumptions on 
the nature of characteristic defects in the formed oxide, which govern the diffusion pro- 
cesses necessary to oxide growth [12, 13]. In turn, this requires a precise knowledge of 
the structure and composition of the ultrathin oxide layers which form spontaneously 



when a bare surface is put in contact with an oxidising environment. These layers 
are often glassy oxides of non-integer stoichiometry, so structural information on the 
oxide lattice or on the defects present cannot be inferred from the corresponding bulk 
oxide structures which form at more elevated temperatures or oxygen pressures [14]. 
Moreover, in the case of ultrathin oxide layers it is very difficult, if not impossible, 
define a precise metal/oxide interface as in the case of thicker crytalline oxides grown 
on metals (see e.g. [15,16]). 

Differently to the growth of oxide layers at high temperature, which is limited 
basically by the thermally activated diffusion of ionic species across the oxide layer, 
the initial reactions of a bare metal surface with an oxidising atmosphere are not 
thermally activated. Namely, spontaneous formation of a few monolayers of oxide has 
been observed at temperatures as low as a few K [17, 18]. Formation of an oxide layer 
takes place if adsorbed oxygen atoms can change place with underlying metal atoms 
and be incorporated below the metal surface. Driving forces for this metal/oxygen 
place-exchange process are in principle available due to the generally high heat of 
adsorption of oxygen on metals, and due to the large electric field which develops at 
the metal/oxide interface as a consequence of electron donation from the metal to 
the incoming oxygen [10, 11]. However, the mechanistic details of the actual exchange 
process are still unknown, and it remains to be understood how this can take place at 
virtually K [6]. 

Our goal is to gain a mechanistic insight into the events of initial oxide formation on 
metal surfaces, prior to the beginning of a diffusion-limited oxide growth process. As 
these events are not thermally activated and occur on a scale of picoseconds to nanosec- 
onds, they are hardly accessible through commonly available experimental techniques. 
Therefore, they will be addressed in this work by means of "first-principles" atomistic 
modelling at the quantum mechanical level [19]. Helped by the fast development of 
accurate and efficient theoretical formalisms and by the ever-increasing power of com- 
puter platforms, this is becoming one of the methods of choice to investigate chemical 
processes in materials science at the atomic scale [20]. In particular, so-called first- 



principles thermodynamics approaches have been widely used in the last few years 
to address the structure and composition of thin oxide layers on metal surfaces (see 
e.g. [14,21,22]). Such approaches give precise information on the different oxide phases 
and structures which are thermodynamically stable at various temperatures and oxy- 
gen pressures. However, low-temperature oxidation reactions may lead to formation 
of intrinsically metastable, often amorphous, oxide structures. Therefore, here we will 
investigate their formation via a first-principles molecular dynamics approach. This 
technique is obviously limited by the small system size and the short simulation time 
addressable, but provide us with a detailed information on the dynamical features of 
the oxidation reactions which would be very difficult to gain otherwise. Moreover, using 
molecular dynamics we are able to obtain structural models of metastable glassy oxide 
layers which may capture the physical and chemical properties of ultrathin oxide layers 
formed at low temperatures more realistically than the corresponding thermodynam- 
ically stable phases. The growth of oxide scales of larger thickness can be simulated 
using atomistic techniques at the classical level, which, however, necessarily need to be 
tuned on previous knowledge from unbiased simulations at the quantum level [23,24]. 

2. Computational details 

All the results presented in the following sections were obtained using molecular dynam- 
ics techniques in which the forces acting on the atoms are computed within the Density- 
Functional Theory (DFT) quantum mechanical formalism. Using a first-principles ap- 
proach, as comprehensively reviewed in [19,20,25,26], is strictly necessary to investigate 
processes involving creation and rupture of covalent bonds. After computing the forces 
within the DFT, the atoms are moved according to the classical equation of motion 
integrated with standard algorithms [27]. All the calculations in the present work were 
carried out in this theoretical framework, using the Lautrec code [28]. Both the min- 
imisation of the electronic states and the dynamics of the atoms were performed using 
the Car-Parrinello (CP) method [29] and the gradient-corrected exchange-correlation 



potential PW91 [30]. Either separable, norm-conserving atomic pseudopotentials [31] 
or Blochl's Projector Augmented Waves approach [32] have been used to describe the 
electron-nucleus interactions including in the simulations only the valence electrons. 
The CP simulations of metallic systems were performed with the algorithm proposed 
in Refs. [33] and [34] . The simulated systems are periodically repeated in all directions 
of space. Surface models were constructed including in the simulation cell a slab of a 
few layers of metal atoms separated from its periodically repeated image by a vacuum 
layer in the direction perpendicular to the surface. The thickness of the slab and of 
the vacuum layer have been chosen so as to avoid spurious effects due to truncation of 
the crystal or interactions with the image systems. Details of the simulation param- 
eters and test calculations to assess the precision of the formalism can be found, in 
particular, in [3,35,36]. 

3. O2 adsorption on bare metal surfaces 

3.1. Experimental background 

Bare metal or semiconductor surfaces are intrinsically very reactive systems. The 
cleavage of a crystal along a given crystallographic direction leaves under-coordinated 
surface atoms and broken, "dangling" bonds which readily participate in chemical 
reactions with molecules in proximity of the surface. A pictorial representation of 
the dangling bonds of a Si(OOl) surface is given in Fig. [T] (top), where the computed 
Electron Localisation Function (ELF) [37] of the system is depicted as a green iso- 
surface. Regions of large ELF far from the atomic cores and of chemical bonds are 
indicative of the presence of lone electrons, and thus of potentially reactive sites for 
chemical reactions. In fact, any Si surface put in contact with the atmosphere at 
normal temperature and pressure conditions would promptly react with dioxygen and 
water molecules and become covered by an amorphous ultrathin hydroxylated oxide 
layer within a few nanoseconds (Fig. [H bottom). Metal surfaces behave in the same 



way, although the electrons available for donation into oxidising molecules are rather 
uniformly delocalised over the atoms composing the surface. 

The reactions between O2 molecules and bare metal surfaces have been studied 
using scanning tunneling microscopy (STM) under ultrahigh vacuum conditions [38- 
40]. When a previously cleaned metal surface is treated with a flow of oxygen gas at 
low pressure, pairs of oxygen atoms separated by distances of about 1 to 3 nearest- 
neighbour metal-metal distances become visible under the STM (Fig. [2]). This suggests 
that O2 reacts dissociatively with metal surfaces by cleavage of the 0-0 bond, leading 
to pairs of single adsorbed oxygen atoms. On the Al(lll) surface at room temperature 
or below, the reaction may proceed with adsorption of only one O atom and ejection 
of the second O atom in the gas-phase [40]. This process is still thermodynamically 
favoured since the heat of adsorption of just a single atom on Al(lll) is more than 7 eV, 
while the energy required to break the 0-0 bond of O2 is slightly larger than 5 eV. 
Both the possibility of O ejection and the large distances (up to ~10 A) which separate 
pairs of atoms after the oxidation reaction at very low temperatures indicate that the 
dissociation process of O2 on the metal surface is associated with a sudden release 
of high kinetic energy. This energy release enables the atoms to overcome relatively 
large diffusion barriers over the surface (of the order of ~0.5 eV) and to cover distances 
much larger than a normal thermally activated diffusion process would allow at the low 
reaction temperatures. Therefore, the term "hot-atom" -mechanism has been used in 
the literature to describe the non-thermally-activated process of dioxygen dissociation 
on metals [38,41]. 

In the following section, first-principles molecular dynamics (FPMD) simulations of 
the process of dioxygen dissociation on several metal surfaces will be reported, focusing 
on the actual driving force for the occurring "hot-atom" dissociation event. 



3.2. "Hot-atom" dissociation of O2 on bare metal surfaces 

Before describing the results of FPMD simulations of the dissociative adsorption of 
dioxygen on metallic surfaces, it is useful to summarise very briefly some features of 
the molecular orbitals of O2 (Fig. [3]). Dioxygen is a paramagnetic molecules, having 
the two degenerate vr* antibonding orbitals each filled with one electron of equal spin. 
The empty a* orbital lies almost 7 eV higher in energy than the vr* orbitals, and 
the Pa; orbital lobes along the direction of the 0-0 bond are evidently compressed 
in the interatomic bond region (see Fig. [31 right). It needs to be noted that within 
our ground state DFT method we are not able to address excited electronic states of 
dioxygen, which may in principles have a role in the oxidation reactions. However, as 
discussed in [35], the mechanisms of onset of oxide formation are expected not to be 
influenced qualitatively by the adiabaticity of the DFT simulations. 

The adsorption of an O2 molecule on the Al(lll) surface has been simulated using 
a periodically repeated surface slab of five (111) atomic layers separated by a vacuum 
layer of the same thickness (Fig. H]). When an oxygen molecule is placed near this 
surface slab, the two unpaired electrons in the tt* molecular orbitals can be visualised 
by imaging of the spin-density of the system. This is defined as the difference Ap 
between the particle densities pa and Pj3 associated with each spin-manifold. In Fig. H^, 
the 0-0 axis of the oxygen molecule is perpendicular to the plane of the page, and 
the radial symmetry of the vr* antibonding orbitals is evident. When the molecule 
is about 3.0 A above the surface, the integrated spin density is about 1.8 electrons, 
indicating a partial charge donation from the surface into the molecule, and thus a 
chemical interaction between molecule and surface. This is also revealed by analysis of 
the local density of states (LDOS) integrated within spheres around the O atoms with 
radii corresponding to the covalent radius of oxygen (Fig. [S]). Namely, the lower edge 
of the LUMO orbital of dioxygen is placed at the Fermi level of the system (Fig. [5l top) 
so that the orbital is not completely empty. As a general shortcoming of standard DFT 
techniques, the barrier for passing from a physisorbed molecular state to a chemically 



adsorbed (and dissociated) state on the Al surface is missing (for a thorough discussion 
of this issue, see e.g. [42]). 

Indeed, after starting the dynamics, the molecule is quickly adsorbed on the surface 
in a process which is characterised by a full quenching of the spin density due to electron 
donation into the half- filled vr* antibonding orbital of O2 (Fig.HJD, c and Fig.E], middle). 
This leads to a gradual increase of the 0-0 distance, until also the a* antibonding 
orbital of the adsorbed molecule becomes partially populated, as revealed by analysis 
of the spin density (Fig. Hid) and of the LDOS (Fig. El bottom). The 0-0 bond 
dissociates as a consequence of this event, the 0-0 distance increasing abruptly in a 
way which is indicative of a non-thermal, "hot-atom" mechanism (Fig. H^, f). The 
sudden release of kinetic energy during the dissociation is large enough to push one 
of the O atoms below the Al surface, while one Al atom is pulled out of the surface 
(Fig. H^). However, in the subsequent dynamics the O atom emerges again above the 
surface and the gained kinetic energy is gradually transferred via the lattice vibrations 
to the whole surface slab. After quenching of the atomic motion, the O atoms are 
stably adsorbed in hollow surface sites, separated by a distance of 5.9 A, i.e. by more 
than two Al-Al distances. 

In a number of similar FPMD simulations, the same "hot-atom" dissociation mech- 
anism has been observed to lead to spontaneous O2 dissociation on Al(lOO), Ti(OOOl), 
Co(OOOl), Cr(llO), TiN(OOl), and notably also on Si(OOl) [36]. We may thus consider 
it as a general mechanism for the initial reaction of dioxygen molecules with bare metal 
or semiconductor surfaces. Among the system investigated so far, only the reaction 
of dioxygen with Pt(lll) resulted in stable adsorption of molecular O2 without dis- 
sociation. In fact, the experimental evidence is that dissociation of O2 on Pt(lll) 
occurs only at temperatures higher than ~150 K [41], while our simulations are all 
performed starting from fully relaxed systems, i.e. at virtually K. However, when 
dissociation is observed experimentally, then the process bears all the typical features 
of a "hot-atom" mechanism [41]. This suggests that thermal motion of the surface is 
necessary to activate the O2 molecule so that the 0-0 distance becomes sufficiently 



large to enable partial filling of the a* antibonding molecular orbital. After that, we 
expect the dissociation to proceed similarly to all other cases considered. 

The dynamical evolution of an O2 molecule during a FPMD simulation of its ad- 
sorption and dissociation on the Al(lOO) surface is reported in Fig. [61 In the top graph 
the evolution of the intramolecular 0-0 distance is plotted together with the evolu- 
tion of the distance between the centre of mass of the two O atoms and the surface. 
The bottom graph reports the evolution of the spin density and the charge on the 
atoms integrated within the Bader regions [43] associated with both O atoms. It is 
interesting to note that after the initial adsorption of the molecule on the surface (left 
vertical dotted line), the charge on the molecule is only slightly more than 1.0 electron. 
Only as a consequence of the "hot-atom" dissociation event (right vertical dotted line), 
which suddenly cleaves the molecular bond, do the O atoms become charged each with 
almost two electrons. This highlights the fact that the hyperthermal dissociation is 
not triggered by Coulomb repulsion between the O atoms, but by the Pauli repulsion 
between the partially occupied p-type orbitals along the 0-0 axis [35]. 

4. Spontaneous formation of ultrathin oxide layers 

4.1. Onset of oxide formation at increasing oxygen coverages 

The same technique used to simulate the adsorption of single oxygen molecules on 
bare surfaces has been applied to investigate further oxidation reactions at increasing 
coverages of adsorbed oxygen on the Al(lll), Ti(OOOl) and TiN(OOl) surfaces. Namely, 
for each system, consecutive FPMD simulations of the adsorption of oxygen molecules 
have been performed starting with a new O2 molecule about 3.0 A above the relaxed 
surface obtained in the previous simulation. On both Al(lll) and Ti(OOOl), O2 is 
again observed to adsorb and dissociate spontaneously with a sudden release of kinetic 
energy. Interestingly, in both systems the process lead to incorporation of O atoms 
underneath the surface layer at a coverage of 0.50 ML. At this coverage, on Al(lll) an 



Al vacancy is created temporarily as a consequence of the "hot-atom" O2 dissociation, 
and is quickly filled by a previously adsorbed O atom nearby (Fig. [71 left). This process 
represents the onset of formation of an oxide phase, as pointed out in Ref. [44]. The 
fact that this is observed in our simulation at a coverage well below the saturation 
of one (1x1) O adlayer is indeed consistent with the experimental finding that oxide 
nucleation can start at coverages as low as 0.2 ML [44]. A similar scenario is observed 
on Ti(OOOl), where the motion of the surface Ti atoms induced by the energy released 
upon dissociation is such that a previously adsorbed O atom is able to penetrate the 
surface layer and bind to the Ti layer underneath (Fig. [TJ middle). These simulations 
strongly suggest that place-exchange processes between metal atoms and O atoms are 
spontaneously activated during the oxidation reaction due to the high kinetic energy 
associated with the cleavage of the 0-0 bond of incoming O2 molecules. 

Interestingly, such place-exchange mechanisms are not observed in simulations of 
oxidation reactions on the TiN(OOl) surface at low oxygen coverages. In this case, 
no O atoms are incorporated underneath the surface, whereas Ti atoms are observed 
to escape from the topmost surface layer (Fig. [TJ right) and bind to the incoming 
oxygen. Also, not all O2 molecules are observed to dissociate spontaneously, indicating 
a decreased reactivity of the partially oxidised surface with respect to the bare TiN 
surface. In particular, increasing the simulation temperature to 600 K was not sufficient 
to dissociate the molecule depicted in Fig. [T]^ right) within 2.5 ps of simulation time. 

4.2. Formation of ultrathin oxide layers 

Further FPMD simulations of reactions with O2 molecules with the partially oxidised 
Al(lll) and TiN(OOl) surfaces were performed to increase the oxygen coverage beyond 
the onset of oxide formation described in the previous section. These are described 
separately for each system in the following two subsections. 
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Formation of amorphous oxide on Al(lll) 

On Al(lll) further reactions with dioxygen proceeded spontaneously at virtually K 
(i.e., both the surface and the molecules were fully relaxed in the initial configuration 
of each FPMD run) up to a coverage of 1.0 ML [35]. After that, O2 molecules initially 
placed at a distance of about 3.0 A above the surface were repelled by the already 
formed ultrathin oxide structure. Therefore, O2 molecules were placed in close contact 
with Al atoms at potentially reactive sites on the outer oxide surface (the initial Al-0 
distances were about 2.0 to 2.5 A). Although some of the molecules were again repelled 
by the surface and desorbed, in a few cases spontaneous dissociation took place, leading 
to growth of the oxide layer up to a coverage of 1.5 ML. The final structure obtained is 
reported in Fig. [Hi Noteworthy is the evident disorder in the oxide structure along the 
direction perpendicular to the surface, while a top- view (bottom left in Fig. [8]) reveals 
a certain order in the surface plane, in which the O atoms are arranged in positions 
roughly correspondent to the hollow sites of the Al(lll) plane. Also interesting is the 
fact that the initially formed oxide structure presents large cavities at the metal/oxide 
interface, which seem to be consistent with the low density of the thin amorphous 
oxide on Al(lll) observed in x-ray photoelectron spectroscopy (XPS) experiments [45]. 
In fact, the XPS evidence suggests that only after reaching a critical thickness does 
the oxide layer undergo a transition to a crystalline state, which is associated with a 
densification of the oxide layer [45]. 

It needs to be stressed that the dynamics of the spontaneous reactions described 
above, and in particular the occurrence of metal/oxygen place-exchange processes after 
dissociative adsorption of O2 molecules, does lead to oxide structures which are not 
necessarily in thermodynamical equilibrium. This is evident by looking at the com- 
puted total energy per unit cell of the oxidised Al(lll) surface at an oxygen coverage 
of 1 ML (the unit cell contains 60 Al atoms and 12 O atoms). The energy of the disor- 
dered structure obtained after consecutive FPMD simulations is ~3 eV higher than the 
energy of an ordered monolayer of adsorbed O atoms on fee surface sites. However, the 
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energy released after the oxidation reactions is very quickly dissipated to vibrational 
modes of the underlying Al lattice, so that the formed disordered structure is not able to 
rearrange back into an ordered ad-layer of oxygen atoms. Moreover, the metal/oxygen 
place-exchange processes are favourable from an electrostatic point of view, since they 
lead to a decrease in the surface dipole which develops as a consequence of donation 
of electronic charge from the metal to the adsorbed oxygen. This can be quantified 
by computing the increase of the work function of the metal surface upon oxidation. 
This increases by 0.32 eV in the case of the spontaneously formed oxide structure and 
by 0.46 eV in the case of the ordered O adlayer. Therefore, a rearrangement of the 
formed disordered structure to an ordered adlayer requires not only thermal activa- 
tion to overcome the diffusion barriers, but an additional activation energy due to the 
unfavourable electrostatics. In conclusion, we expect such rearrangements not to take 
place at room temperature or below, while they are certainly not excluded at higher 
oxidation temperatures (at which, however, other ordered ultrathin oxide structures 
may become thermodynamically favourable with respect to an O adlayer [46]). 

Formation of ordered oxide on TiN(OOl) 

The final snapshots of six consecutive FPMD simulations of the initial oxidation of a 
TiN(OOl) surface model containing 9 Ti and 9 N atoms in the periodically repeated 
simulation surface cell are reported in Fig. O As outlined above, the initial reactions of 
O2 molecules with the TiN(OOl) surface take place without incorporation of O atoms 
underneath the surface. Instead, Ti atoms are observed to escape the surface plane to 
become coordinated by the incoming oxygen atoms, leaving behind Ti vacancies in the 
TiN lattice (Fig. [Hb^e). In this process, the surface Ti atoms leave their coordination 
shell of five nitrogen atoms as a consequence of their exceedingly high formal oxidation 
state when oxygen atoms are adsorbed on the surface [3]. The process can also be 
understood in terms of an increase of the ionic character of the Ti-N bonds leading to 
their destabilisation after the formation of the Ti-0 bonds [47]. Therefore, the initial 
oxide formation on TiN appears to be driven by a chemical driving force which induce 
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the formation of a mixed- valence TiN^-Oy surface phase [3]. 

In many of the simulations, dioxygen dissociation did not occur spontaneously, so 
that O2 molecules remained stably adsorbed on the partially oxidised surface for the 
entirety of the simulation runs. In the last simulation, the temperature was gradually 
increased to about 600 K within the first 2 ps of simulated time, then the system 
was annealed at this temperature for about 2 ps. Subsequently the temperature was 
gradually increased to about 900 K within 1 ps and the system annealed for about 
4 ps before quenching the atomic motion and fully minimising the atomic position. An 
interesting feature observed during this simulation is the diffusion of Ti vacancies to 
the second atomic layer below the surface and the creation of Ti interstitials between 
the first and second Ti layers and at the metal/oxide interface. Moreover, the vacancies 
created seem to have a tendency to cluster together, as can be seen in Fig. [TD] where 
the final structure obtained is reported showing two periodically repeated simulation 
cells. Clustering of vacancies at the metal/oxide interface has been indeed observed 
experimentally upon selective Al oxidation of an intermetallic TiAl surface [48] . 

During the dynamics, the initially present O2 molecules dissociated and the oxide 
structure underwent a profound rearrangement to form an infinite chain of fourfold- 
coordinated Ti atoms (Fig. [TOl) above the existing oxide layer. This structure resem- 
bles very closely the proposed "ad-molecule" reconstruction of anatase TiO2(001) [49]. 
Other Ti atoms remained in a fivefold coordination (being bound both to oxygen atoms 
and to nitrogen atoms of the substrate), as commonly found on titanium dioxide sur- 
faces, such as, e.g. rutile TiO2(110) [50]. The ordering achieved in the final structure is 
remarkable given the short simulation time accessible to our technique. It is probably 
due to the strict matching conditions imposed by the underlying TiN lattice, on which 
O atoms can only bound to Ti sites, and Ti atoms only to N sites. Notably, at this 
stage it is already possible to identify three chemically different Ti species. Namely, 
(i) the fourfold-coordinated Ti atoms bound to only oxygen, (ii) the Ti atoms at the 
interface bound to both oxygen and nitrogen, and (iii) the Ti atoms of the bulk bound 
to only nitrogen atoms. These three chemical species are indeed typically found in 
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XPS investigations of the oxide layers grown on TiN substrates [3,47,51]. 

Finally, it is to be noted that in spite of the profound rearrangements of the oxide 
layer formed, the N atoms of the TiN lattice remained in their equilibrium positions. 
Therefore, formation of thicker Ti02 layer and release of gaseous N2 are expected to 
involve a slower, activated process. This may be consistent with a number of XPS 
experiments where TiN^O^ compounds have been observed to form uniformly over the 
entire region of a TiN film in the early oxidation stages, prior to the formation of an 
amorphous Ti02 layer. This is then converted into a crystalline Ti02 film only after 
a thermal treatment at higher temperatures [47,51]. Moreover, the native oxide layer 
formed at room temperature has been found to be substantially thinner on TiN than 
on Ti [51], which reflects the high stability of the TiN^-Oj, surface phase. 

5. Conclusions and outlook 

Despite evident differences among the systems investigated, all our FPMD simulations 
indicate that the onset of oxide formation on metal surfaces at low temperature is driven 
by the kinetic energy released during the dissociative adsorption of O2 molecules. Af- 
ter the single reactions, the oxide structures obtained become quickly frozen into a 
metastable state due to energy dissipation via lattice vibration before the next molec- 
ular reaction can take place. In fact, the rate of collision r of a molecular gas with a 
solid surface is: 

r= , ^ (1) 



/27rMA;BT 

where M is the molecular mass, k-^ the Boltzmann constant, T the absolute temperature 
and p the gas pressure. From this is follows that at a pressure of 1 atmosphere and a 
temperature of 300 K, 1 nm^ of metal surface would be hit by about 3 oxygen molecules 
every nanosecond. This justifies our choice of quenching the atomic motion at the end 
of each simulation (lasting no more than 10 ps) and starting the next consecutive 
simulation of O2 adsorption from fully minimised atomic positions. 
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Provided that all molecular collisions lead to dissociative oxygen adsorption, as 
can be indeed observed on many metal surfaces, this would lead to formation of a 
monolayer of oxide within 2 to 3 ns. Only after formation of a few ML of oxide does 
the probability of reaction with further oxygen molecules rapidly decrease, in a way 
which is strongly dependent on the metal considered [8]. The case of Al(lll) is a 
peculiar one, since the fraction of oxygen molecules which actually react upon collision 
(the so-called sticking coefficient) is relatively low (10^^ to 10^^) at very low oxygen 
coverages [52]. Therefore, an ultrathin oxide layer on Al(lll) is expected to form in 
less than a microsecond under normal atmospheric conditions. 

The emerging scenario is one in which surfaces exposed to an oxidising atmosphere 
become almost instantaneously covered by an ultrathin oxide layer whose structure and 
composition may not represent the thermodynamical equilibrium situation. However, 
this does not mean that amorphous thin oxide layers are necessarily metastable with 
respect to the corresponding crystalline phases. As outlined in [13,53], factors such as 
lattice mismatch, stress in the oxide layer, formation of dislocations at the metal/oxide 
interface etc. will have a strong influence on the actual stability of the oxides formed. 
At the level of the simulations described in this work, the formation of vacancies and 
their clustering underneath the oxide layer may have a direct influence on the stability 
and on the physical properties of the native oxide layers. This is an issue of particular 
importance during the formation of native oxides on alloys via the selective oxidation 
of only one element present in the alloy, and will be the subject of future investigations. 

Another factor that should be taken into account is the presence of water vapour in 
the atmosphere at normal conditions, which may condensate as a thin film on oxidised 
metal surfaces. Water has been observed indeed to readily react dissociatively with an 
oxidised Si(OOl) surface in FPMD simulations [54]. Thus, in general the oxide layer 
will be terminated by an amount of hydroxyl groups, similarly to what is shown in the 
bottom panel of Fig. [H 

Our future efforts will be directed at obtaining models of native oxide layers on tech- 
nological metal surfaces capturing their chemical and physical properties under normal 
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atmospheric conditions in a possibly realistic way. Knowing the essential details of the 
structure, composition and electronic structure of native oxide layers may then enable 
us to construct analytic potentials to study the interactions between odixised surfaces 
and the external environment at the classical simulation level, or within hybrid quan- 
tum/classical formalisms. This will certainly be necessary in order to address processes 
which occur on larger size and time scales than what is within the reach of purely quan- 
tum mechanical schemes. A prototypical example would be the initial adsorption on 
metal surfaces of proteins which mediate the interactions between biomedical implants 
and living cells, just to name a field which is of particular interest both to the author 
and to whom the present journal issue is dedicated. 
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Figure 1: Top: Model of a bare reconstructed Si(OOl) surface displayed with its com- 
puted Electron Localisation Function (green isosurface) showing the presence of reac- 
tive dangling bonds on exposed Si atoms. Bottom: model of an oxidised and hydrox- 
ylated Si(OOl) surface as obtained in a series of FPMD simulations [36,54]. 
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Figure 2: Left: STM image of a bare Al(lll) surface upon exposure to an oxygen gas 
at 150 K. Distances between pairs of oxygen atoms are indicated in units of the Al-Al 
nearest neighbour distance. Right: Histogram showing the relative number of atom 
pairs separated by the Al-Al distances reported on the abscissa. The dark spots and 
the histogram line indicated with "1" are in fact associated with isolated oxygen atoms 
produced via non-adiabatic abstractive adsorption processes [40]. (Reproduced with 
permission from Ref. [39], copyright Elsevier 2001). 
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Figure 3: Scheme of the valence molecular orbitals of O2 (left) in its ground state and 
their corresponding associated particle densities (right) represented as blue and green 
isosurfaces for filled and empty orbitals, respectively. 
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Figure 4: Snapshots from a FPMD simulation of the adsorption of a dioxygen molecule 
on the Al(lll) surface [35]. The spin-density of the system is represented as an orange 
isosurface. Note the partial filling of the a* antibonding orbital (snapshot d) which 
triggers the "hot-atom" dissociation event. 
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Figure 5: Evolution of the computed local density of states (LDOS) projected onto 
the O atoms during the adsorption of an O2 molecule on the Al(lll) surface. The 
continuous and dashed lines correspond to the majority and minority spin manifolds, 
respectively, and the Fermi level is indicated with a vertical line. From top to bottom, 
the graphs refer to an oxygen molecule far from the surface, adsorbed on the surface 
(filling of the vr* orbitals and spin-quenching), and at the initial moment of dissociation 
(partial filling of the a* orbital) . 
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Figure 6: Dynamics of an O2 molecule adsorbing dissociatively on the Al(lOO) surface 
simulated by FPMD. Top graph: evolution of the 0-0 distance and of the height of 
the centre of mass of the O atoms with respect to the average height of the surface Al 
atoms. Bottom graph: evolution of the spin-density and of the atomic charge integrated 
within the Bader regions associated with both O atoms. 
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Figure 7: Structures of partially oxidised metal surfaces as obtained in FPMD simu- 
lations (only the three topmost metal layers are shown, for clarity). The oxygen atom 
coverage of the surfaces is 0.5 in the case of Al(lll) and Ti(lll) and 0.67 in the case 
of TiN(OOl). Note the spontaneous incorporation of O atoms underneath the surface 
in the case of Al and Ti, the extraction of metal atoms from the surface layer in the 
case of Al and TiN, and the presence of an adsorbed O2 molecule on TiN. 
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Figure 8: Model of a native oxide layer grown on the Al(lll) surface (the oxygen cov- 
erage is 1.5 ML), as obtained in a series of FPMD simulations. Orthogonal projections 
and perspective view (top right). 
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Figure 9: Final snapshots of consecutive FPMD simulations of the initial oxide forma- 
tion on the TiN(OOl) surface. The oxygen covrage increases from (a) 0.22 ML to (f) 
1.33 ML. 
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Figure 10: Model of a native oxide layer grown on the TiN(OOl) surface (the oxygen cov- 
erage is 1.33 ML), as obtained in a series of FPMD simulations. A fourfold-coordinated 
Ti atom is indicated with Ti(4), a fivefold-coordinated Ti atom with Ti(5), and an atom 
in interstitial position with respect to the TiN lattice with Ti(I). 
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